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1. Introduction: theory and computational setup 

In the last few years, the ETM Collaboration (ETMC) has been carrying out state-of-the-art 
lattice QCD simulations with Nf = 2 dynamical flavours (sea quarks) and 'lightish" pseudoscalar 
meson masses (300 MeV < nips < 550 MeV). Strangeness clearly enters the game in a partially 
quenched context. Several physical quantities are currently being analyzed for a few lattice spac- 
ings. In the present work we will present preliminary results on: (i) the vector meson mass mK* , 
(ii) its decay constant /k*, (iii) the ratio of tensor to vector decay constants fr/fv, computed at the 
K* physical mass and (iv) the Bk parameter for neutral Kaon oscillations. Other subgroups of the 
collaboration are working on decay constants in the light and strange quark sector [JlJ, ^]. 

ETMC simulations are performed with the tree-level Symanzik improved gauge action. The 
Nf = 2 sea quark flavours are regularized by the standard Wilson fermion action with a maximally 
twisted mass term (referred to as "standard tmQCD case") . This means that the the two light 
flavours are organized in a flavour doublet iff^ = (m d) and the fermion lattice action is given by 



■^tm — W 



Dw + il^qT^Ys 



¥ (1-1) 



with Dw denoting the critical Wilson-Dirac operator. This formulation has well known advan- 
tages, amongst which: (i) Renormalization properties are much simpler than with standard Wilson 
quarks, in many cases of interest (e.g. pseudoscalar decay constants, chiral condensate, 8^...); (ii) 
improvement is automatic at maximal twist [Q]. One should note, however that for Weak Matrix 
Elements (WMEs) of 4-fermion operators (e.g. B^), it is not possible to retain both automatic 
improvement (through the standard tmQCD formalism, with all flavours at maximal twist) and 
multiplicative renormalization for the relevant operator; see refs. [Q, |6|]. One way out is provided 
by the Osterwalder-Seiler (OS) variant of tmQCD, in which valence quarks enter with a distinct 
action for each flavour and a fully twisted mass term [0] : 

= Vf Dw + i^fYs ¥f f = u,d,s--- (1.2) 

Suitable combinations of Hf signs for each flavour in the above action ensure automatic improve- 
ment and multiplicative renormalization for say, Bk- The OS option is a compromise, since uni- 
tarity issues arise at finite lattice spacing when sea and valence flavours are treated differently. 
Nevertheless, in our partially quenched setup (Nf = 2 sea quark flavours and a valence strange 
quark) compromising unitarity is anyway unavoidable in any regularization. 

The ETMC runs are performed at three gauge couplings j3. Here we report work in progress, 
confining ourselves to the "master run" at j3 = 3.90, corresponding to a lattice spacing of a ~ 
0.086(l)fm (i.e. a^^ ~ 2.3GeV). The lattice volume isV = 24^ x 48. Our ensemble consists of 
240 gauge field configurations for the /T* -meson mass and decay constants and 200 configurations 
for Bk- The sea quark mass is set at five values = 0.0040, 0.0064, 0.0085, 0.0100, 0.0150, 
corresponding to pseudoscalar meson masses in the range 300MeV < mps < 550MeV. The valence 
quark masses are set at seven values; those of the sea quarks plus the values ajJ. = 0.0220,0.0270, 
which are meant to facilitate the interpolation to the physical strange quark mass. We use the 
calibration results of previous ETMC work [Q, ^; e.g. the physical light quark mass is at a^d = 
aiJL{mji) = 0.00079, whereas the physical strange quark mass is at apLs = aiJ.{m^) = 0.0217(10). 
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For Sk only, we checked for finite volume effects by running also at V = 32^^ x 64 at the lowest 
quark mass (a/i = 0.0040). 

For the computation of correlation functions, we used stochastic sources of the extended "one- 
end trick" of refs. [|Io|,[ri[]. For a concise exposition of the method see also ref. [12]. 



2. Vector meson masses and decay constants 

The vector and tensor decay constants are defined by the formulae: 

{0\%\V;X)=fvetmv {Q\%k\V;X) = -ihe^mv (2.1) 

Vector meson masses and the decay constants /y, fj are computed from the correlators (in contin- 
uum notation) Cy,i/,(0 = ^x,k < M^) > and CTo,To,{t) = Lri < ^kix) ^\k{0) >. The 
best estimate of the vector meson effective mass is provided by the correlator Cy^Vk (0- The estimate 
from the correlator Cj^T^it) is much noisier. We compute f'^ / f^\K* from the large time asymptotic 
behaviour of the ratio [CT^^T^tit) /CvtVk{t)Y'^ ■ 

As discussed in ref. [jl3|], the correctly normalized vector current in the tmQCD fermion regu- 
larization is 1^ = Z^AjJ", while in the OS setup we have "V^ = ZyV^^ . With f^y = s^vpaTpa, we 
similarly have for the tensor density = ZtT^^ and = ^tT^^v- We use the RI/MOM 



estimates for the (re)normalisation constants, computed in a tmQCD framework in ref. [13], ob- 
taining Zr(a-^i3 = 3.9) = 0.769(4), ZA(j8 = 3.9) = 0.771(4) and Zy(j8 =3.9) =0.6104(2). The 



slight differences between the above values and the ones of ref. [13] are due to a revised analysis. 

We wish to highlight straightaway the two problems we have encountered: (i) For all sea quark 
masses, when the valence quark masses are in the lightest range (say, around /^vai = 0.0040), the 
vector meson effective mass has a poor plateau. The situation already improves at the next mass 
value of our simulation, /^vai = 0.0064. Nevertheless, since the signal-to-noise ratio behaves as 
expected (i.e. it drops like exp[— (my — nips)f]), the p-meson mass and decay constant can still be 
extracted (see results presented in ref. (ii) A poor quality vector meson effective mass is also 
seen when pLvai < }lsea- This problem is absent in the pseudoscalar channel. 

We thus proceed as follows: the strange vector meson consists of two valence quarks (/i/,/i/i). 
We compute the necessary observables (vector meson mass, vector decay constant and the ratio of 
the tensor to vector decay constant) for all combinations of apLi = apLgea and a/X/, = 0.0150, 0.0220, 
0.0270. In this way unitarity holds in the light quark sector, while the heavy valence quark mass, in 
a partially quenched rationale, spans a range around the physical value /i,. Examples of the quality 
of our signal (for the tightest a/i/ = ajXsea masses) are given in Fig. |l[ 

Subsequently, keeping ajih fixed, we extrapolate these results linearly to the physical down 
quark value apL^ = 0.00079 (alternative fits are currently under investigation). This is repeated for 
all three a/i/j values, before the results are linearly interpolated to the physical strange quark mass 
apis = 0.0217(10). Examples of the quality of such extrapolations and interpolations are displayed 
in Fig. ^ Our final results are (tmQCD regularization on the Ihs and OS regularization on the rhs) 

aMv\K* =0.422(10)(03) aMv\K* = 0.437(08) (04) 
afv\K* =0.106(05) (02) afv\K* = 0.117(03)(01) (2.2) 
h/MK' =0.770(20) (09) h/MK' =0.759(19)(07) 
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Figure 1: Vector meson quantities plotted as a function of the operators' time separation, computed at 
a^i = fljUieu = 0.0040 and a^U/, — 0.0150. (a) effective mass; (b) decay constant; (c) ratio /r/ fy. The quoted 
results are obtained by taking data over the time plateau 11 < xq < 16. 



The first error is statistical. The second is systematic and includes the uncertainty of the apis 
estimate plus (for the decay constants) that of the renormalization factors Zv,Za,Zj. The fairly 
good agreement between tmQCD and OS estimates suggests that cutoff effects in the vector channel 



are small, in accordance with the expectations of ref. |15]. This result must clearly be confirmed 
by repeating the analysis at other j8 values. The experimental values of the vector meson mass and 
decay constants, expressed in lattice units (at jS = 3.9, we use a ~ 0.086 fm) are artiK* = 0.381 and 
fl/if* = 0.0927, which are reasonably close to our lattice estimates. Finally, the tensor to vector 
decay constant ratio is RG-run from the scale a^' ~ 2.3GeV to the usual reference scale of 2GeV, 
obtaining fr / fvW* = 0.764(19)(03). This is pretty close to the continuum limit quenched result 
/r/Ak' =0.74(2) of ref. 

3. The K-meson bag parameter 

We now apply the proposal of ref. in an O(fl') -improved calculation of Bk- Recall that, 
since we require both automatic improvement and multiplicative renormalization of this quantity, 
the setup is "mixed": the two Kaon valence quarks are maximally twisted in standard tmQCD 
fashion, while the two anti-Kaon ones are OS valence quarks (or vice versa). In this way, the 
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Figure 2: (a) Extrapolation of the vector meson decay constant fr to the physical down quark mass an^ , 
at fixed heavy quark mass a/X/, = 0.0150. (b) Interpolation of fr{apLd) to the physical strange quark mass 
alls of the vector meson decay constant fr. The red circles correspond to data points, the blue squares to the 
extrapolation and interpolation results. 



continuum (renormalized) four-fermion operator is related to the lattice bare one as follows (the 
notation should be obvious) 



(3.1) 



With xq = w m arbitrary reference time-slice, the relevant 3-point correlation function consists 
of two "K-meson walls" with noise sources at fixed times (w and w + T /2) and a 4-fermion operator 
living at different time-slices xq, with w < xq < w + r/2. The signal quality is improved in various 
ways. For example, the reference time w is varied from configuration to configuration. Moreover, 
as already stated previously, the valence quarks emanating from w are of the OS variety, while 
those emanating from w + T /2 are standard tmQCD ones. This situation is also reversed and the 
corresponding correlation functions are suitably averaged. 

The 5K-parameter is measured from correlation functions with the light valence quark mass 
kept equal to the sea quark mass {a}Xi = a/isea)- At fixed heavy valence quark mass ajXh, we fit 
the light mass behaviour in ajXi, using the SU(2) Partially Quenched Chiral Perturbation Theory 



(PQ-;t;PT) result of refs. [Il7|,[18|]: 

(3.2) 



B{Hh) = B^{ni,) 1 + b{Hh)—Hi - ^^^^ ^^ ixi In f 



f2 ^' 32712/2 ^^-"V ^2 

With the coefficient IBq//^ known from earlier chiral fits of the light quark sector ^ ^], the above 
relation requires a two-parameter fit (B^^ and b{pih)) in the chiral region. In this way we extrapolate 
B{iih) at the physical down quark mass apLd- The result is shown in Fig. ^a). Alternative polyno- 
mial fits are currently under study. The final step is to linearly interpolate the estimates to 
the physical strange quark mass value ajis- This is shown in Fig. ^(b). Note that the result of the 
highest ajXh has not been included, as it lies rather far from apLg- 

Finite size effects appear to be under control, since at pLsea = 0.0040 and Hi = Hf, = 0.0100, 
we find B^^^^ = 0.591 (5) for the L = 24 lattice and B^f?'"" = 0.598(8) for the L = 32 one. 
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(a) (b) 

Figure 3: (a) Extrapolation of Z?k to the physical down quark mass AjUj, at fixed heavy quark mass a/i/, = 
0.0220. The Bk datapoints are computed at degenerate light and sea quark masses, jli = /isea and fitted to 
eq. (3.2). (b) Linear interpolation of Bk, as a function of the heavy valence quark mass a/i/,, to the physical 
strange quark mass a^s- 



Renormalization is again performed in the RI/MOM scheme [19]. The quality of our prelimi- 
nary results for the multiplicative renormalization factor Zva+av is shown in Fig. ^(a). The absense 
of mixing with "wrong chirality" operators is explicitly demonstrated in Fig. ^^b), where all mixing 
coefficients are shown to be vanishing. We estimate Zva+av{^ = 3.9;2GeV) = 0.454(18). 
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Figure 4: (a) RI/MOM computation of the multiplicative renormalization factor Zva+av- (b) Mixing coeffi- 
cients Ai^ik— 1 , • • • ,4) with other four-fermion operators with "wrong chirality". 



The results reported in this work are very encouraging, but a word of caution is in place here. 
At fixed j3, the two Kaon states, obtained with different regularizations (i.e. tmQCD and OS) are 
not degenerate, differing by 0{a^) discretization terms. The two different exponential decays, as 
well as the factors m'^ and of the matrix element < K'^{m}-^)\OvA+Av\K-^{ni^) m}^m^, 
cancel out in the ratio of Bk- The Kaon mass splitting is sometimes quite significant; for the 
j3 = 3.90 case in hand, with a/i/ = a/isea = 0.0040 and ajih = 0.0220 we have m'jf = 0.2391(07) 
and = 0.2923(16). It is important to monitor the size of this splitting with increasing /3, in 
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order to confirm that it vanishes Uke a^. 

Our preliminary result at a single lattice spacing is BK(2GeV;RI/MOM) = 0.56(2), corre- 
sponding to = 0.77(3). 
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